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(57) Abstract: A method and apparatus of generating a hybrid three dimensional reconstruction of a vascular structure affected by 
periodic motion comprises placing an object (50) affected by periodic motion to be imaged in an imaging region of an x-ray system 
22, the object having a vascular structure. At least two x-ray images of the vascular structure are acquired (104, 204). Indicia of the 
phases of periodic motion are obtained (104, 52) and are correlated with each of the x-ray images. At least two x-ray images from 
a similar phase of periodic motion are selected (108). A three dimensional modeled segment of a region of interest in the vascular 
~ structure is generated (1 10, 210), the modeled segment reconstructed using the selected x-ray images from a similar phase of periodic 
motion and the region of interest only a portion of the imaged vascular structure. A three dimensional volumetric reconstruction of 
a vascular structure is generated (1 12, 212, 207) that is larger than the modeled segment. The modeled segment of interest (148) and 
}^ the volumetric reconstruction of the larger vascular structure are combined and displayed (220) in human readable form. 
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Rotational Angiography Based Hybrid 3-D 
Reconstruction of Coronary Arterial Structure 



The present invention relates to x-ray coronary 
5 angiography and is particularly related to an apparatus that 
applies rotational angiography for acquiring views to 
generate a hybrid reconstructed and modeled representation. 
The present invention finds particular application in 
conjunction with diagnostic medical imaging for use in 
10 cardiac catheterization, diagnosis and interventional 

treatment and will be described with particular respect 
thereto . 

Coronary artery disease remains a major cause of 
morbidity and mortality in the United States. In an effort to 

15 reduce the mortality rate, there has been a marked increase 
in the number of catheterizations procedures performed. In 
the future, the number of catheterizations is expected to 
increase due to recent advances in stent technology (e.g. 
eluting stent) and imaging capabilities . 

20 Coronary catheterizations are generally performed using 

X-ray angiography. Choosing the correct stent dimensions for 
use during the interventional procedure is often difficult 
using the traditional 2D projection images due to vessel 
foreshortening and overlap. 

25 Based on the specific anatomy of a specific coronary 

branch of clinical interest and skill of the clinician, 
several routine 2D single plane or biplane angiographic 
images from different arbitrary viewing angles are acquired 
in order to provide a working image for the clinician to 

30 visually estimate and derive the length and diameter of the 
stent to be used for subsequent intervention. Once the 
clinician decides upon the stent dimensions, a working view 
is chosen to position the x-ray system for the clinician to 
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use during the interventional procedure. In this working 
view, it is desired that the vessel segment of interest is 
the least foreshortened and is not blocked from sight blocked 
from view by overlapping structures e.g. other coronary 
5 arteries, implants or patient structures- If the working view 
selected based on the experience of the clinician is not as 
useful as desired, an alternate view is selected. 

Since the patient receives a dose of contrast agent or 
medium and radiation for each view the clinician chooses to 
10 evaluate, this trial-and-error method of selection of 

appropriate viewing angles potentially exposes the patient to 
large amount of contrast medium (dye) and radiation, even 
before the actual interventional treatment of the patient 
begins*. . The contrast media presently in use places stress on 
15 a patient's kidneys and it is desirable to limit the dose 
amount of contrast medium to a patient. In patients that 
have reduced kidney function or disease, the reduction of 
contrast medium dose is even more important. In addition, 
the contrast media can harm the heart and is limited to 
20 prevent permanent damage- Likewise, it is desirable to 

reduce the patient's exposure to x-rays during the procedure. 

In order to improve the assessment and treatment for 
patients relative to the ^trial-and-error' approach, several 
researchers have proposed various methods to construct a 
25 three-dimensional surface model of the entire arterial trees 
for the coronaries using two or three static 2D x-ray image 
acquisitions from different arbitrary viewing angles, 
generally needing considerable clinician experience, judgment 
and interaction. As a result, after obtaining the static 

30 acquisitions, the clinician has the ability to use the 2D 
views to generate and view a model of the entire three- 
dimensional coronary tree from any angle- By manipulating 
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the model in the image the clinician is able to choose an 
^optimal working view' without the use of extra radiation or 
dye. Modeling the entire 3-D coronary tree is still time 
consuming and generation of the coronary artery tree 
5 reconstruction and model are deleteriously affected by the 

periodic cardiac motion. In addition the few selected views 
for model construction are limited to those acquired based on 
the skill of the clinician without the benefit of knowing the 
unique structure for a particular patient- The individual 

10 variation in coronary tree structure may decrease the value 
of one of the standard static views thereby requiring 
acquisition of yet another image and corresponding dose of 
contrast medium. 

In the last few years, rotational angiography (RA) has 

15 proven to be a very accurate and effective diagnostic tool in 
the treatment of cerebral vessel malformations. In this 
approach, the C-arm rotates rapidly around the patient's head 
while several X-ray projections are acquired. The 
reconstructed cerebral vessels can be viewed from different 

20 viewing angles, while only one contrast injection is given . 
Due to the high reproducibility of the rotational 
acquisitions, the fast rotation speed, and the static nature 
of the cerebral vessels, the projections can be used for 
volumetric reconstruction providing very high detail and 

25 accuracy. In addition, the process is automated for the 
stationary cerebral vessel structures . 

It has been demonstrated that the use of rotational 
arigiography for coronary vessel acquisitions yields better 
stenosis severity estimations and reveals lesions that were 

30 missed by only applying a few of the traditional static image 
acquisitions. However, straightforward volumetric 
reconstruction of coronary arteries has several difficulties. 
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First, due to the beating of the heart and patient 

respiratory motion, present techniques only yield a rough 

representation of the coronaries. Given the need for accurate 

information regarding the dimensions of the coronary lesion 

5 for selection of the proper stent or other interventional 

device or procedure, it is desirable, to have more precise 

measurements of the coronary arterial tree than possible with 
these rough representations . More precise information 

improves the diagnostic information available to clinicians 

10 in selection the optimum treatment for the patient. Although 
recent developments with alternative reconstruction schemes 
are very promising, an accurate reconstruction that can 
directly be used for the dimensioning of the stenoses in • 
hioman subjects is not yet available or reliable enough for 

15 clinical use. Even if a suitable volumetric reconstruction 

could be performed, additional user interaction is needed to 
arrive at the correct stent dimensions. 

Some efforts have been made to improve the rough 
reconstructions and images by attempting to compensate the 

20 rotational angiograph for periodic cardiac motion by gating 
the acquisition of the images to a predetermined portion of 
the cardiac cycle to provide a more desirable image for the 
reconstruction. For example. International Patent 
Application No, WO 02/3 6011 entitled ^^Method and Apparatus 

25 for 3D-Rotational X-Ray Imaging" describes using an x-ray 

imaging method for a 3D reconstruction of a coronary vascular 
system with a relatively long run of a scan rotation over 15 
sec. to preferably 20 sec. and limited to covering 10° per 
second. The images are acquired using ECG cardiac triggering 

30 so that the images are obtained at the same point in the 
cardiac cycle. 
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In U.S. Patent No. 6,324,254 entitled ^^Method and X-ray 
Device for Picking Up X-Ray Images of a Substantially 
Rhythmically Moving Vessel or Organ" describes an x-ray 
device that is slowly moved along a circular orbit at an 
5 angular velocity of less than 6° per second, during which a 
nimber of digital x-ray images are picked up, with the image 
pick up being triggered by the vessel motion or organ motion, 
e.g. an EGG device. The device rotates through 150°- 200° 
which indicates a scan time of 25 sec. to 33 sec. The patent 
10 describes using even slower rotation rates of 2° and 0.5° 
which indicate a potential total scan time between 75 sec. 
and 400 sec. The increased scan time is necessary to obtain a 
sufficient number of data acquisitions for the same phase of 
heart cycle. 

15 During a typical 6 sec. C-arm scan, only 30 images are 

obtained at the same phase of the heart cycle. Using this 
low number of images for the reconstruction results in poor 
quality of the reconstructed image. It has been indicated in 
the above reference that 100 such images at the same phase of 

20 the heart cycle are desired to reduce noise in the resulting 
reconstructed 3-D volume - 

In addition, even with the improved reconstructed images 
using the required longer scan time, the resolution of the 
arterial trees still leaves room for improvement in 

25 determining the Quantitative Coronary Analysis for the 

dimensions of a stenosis or region of interest for treatment 
by the clinician. Furthermore, additional time consuming 
clinical user interaction is necessary in order to extract 
the correct dimensions of the lesion from the volumetric 

30 reconstruction. 

However, due to the teaching in these disclosures 
requiring lengthening the rotational scan time, both of these 
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techniques exacerbate the substantial problems of increasing 
patient dose for both of the contrast medium and x-ray 
exposure. As a result some patients with existing kidney- 
problems may not be able to withstand the additional required 
5 contrast medium dose for an improved reconstructed view of 
the coronary tree. This will reduce the quality of 
information available to the clinician in treatment of these 
patients. In addition, all patients would receive larger 
doses of x-radiation during these extended scan times. 

10 Therefore it is desirable to have an angiographic 

procedure that will result in generating an improved hybrid 
three dimensional coronary arterial tree compensated for 
cardiac motion that provides (i) accurate modeled three 
dimensional representations of selected segments of the 

15 coronary arterial structure suitable for Quantitative 

Coronary Analysis, (ii) a hybrid modeled and volumetric 
reconstruction of the coronary arterial tree for selecting an 
improved optimal working view. It is also desirable to 
provide an accurate modeled representation of selected 

20 segments of the arterial tree to reduce foreshortening and 
overlap so that an optimal view of the lesion may be 
determined in a clinically practical time period for 
diagnostic and treatment use by the clinician during the 
interventional aspect of the patient's treatment. 

25 The present invention is directed to a method and 

apparatus that satisfies the need to provide a clinically 
useful optimal view of the coronary arterial tree for 
coronary interventional procedures while improving the 
clinicians ability to obtain accurate measurements of a 

30 stenosis or lesion and still reducing the amount of time to 
generate clinically useful reconstructed and modeled 
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segments, contrast agent and radiation dose to which the 
patient is exposed. 

A method is presented to assist the clinician in plan- 
ning a diagnostic or interventional procedure while the 
5 patient is already on the catheterization table. Based on 
several selected projections from a rotational X-ray 
acquisition, both a volumetric cone-beam reconstruction of 
the coronary tree as well as a three-dimensional surface 
model of the selected vessel segment of interest is 

10 generated, i.e. a partial model of the coronary arterial tree 
with the segment of interest. The proposed method provides 
the clinician with clinically accurate length and diameters 
of the vessel segment of interest as well as with an optimal 
working view. In this view, the gantry is positioned such 

15 that, for the vessel segment of interest, the foreshortening 
and vessel overlap are reduced during the entire heart cycle. 

A method in accordance with principles of the present 
invention includes a method of generating a hybrid three 
dimensional reconstruction of a vascular structure affected 

20 by periodic motion. The method comprises placing an object 

affected by periodic motion to be imaged in an imaging region 
of an x-ray system, the object having a vascular structure 
and acquiring at least two x-ray images of the vascular 
structure. Indicia are obtained of the phases of periodic 

25 motion and correlating the indicia with each of the at least 
two x-ray images. At least two x-ray images are selected 
from a similar phase of periodic motion and a three 
dimensional modeled segment of a region of interest in the 
vascular structure is generated. The modeled segment is 

30 reconstructed using the selected x-ray images from a similar 
phase of periodic motion and the region of interest is only a 
portion of the imaged vascular structure. A three 
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dimensional voliimetric reconstruction of a vascular structure 
larger than the modeled segment is generated and is combined - 
with the modeled segment of interest. The combined 
reconstructed vascular model and volumetric reconstruction is 
5 displayed in human readable form. 

Another more limited aspect of the method applying 
principles of the present invention includes using x-ray 
images from all of the phases of the periodic motion for 
generating the volumetric reconstruction includes . In 

10 another limited aspect in accordance with the present 

invention the volumetric reconstruction is generated with 
gated images from a similar phase of periodic motion. 

In another limited aspect of the present invention, the 
volumetric reconstruction includes a first portion of the 

15 vascular structure reconstructed using x-ray images from all 
of the phases of periodic motion and a second portion of the 
vascular structure reconstructed with gated images from a 
similar phase of the periodic motion. A more limited aspect 
of the method includes combining the reconstruction of both 

20 of the first portion and second portion of the vascular 
structure with the modeled segment of interest. 

Another limited aspect in accordance with principles of 
the present invention is that the obtained indicia of phases 
of periodic motion is representative of cardiac motion and is 

25 provided by an ECG signal. 

Another limited aspect in accordance with the present 
invention includes computation of an overlap map, the overlap 
map is computed by integrating all the gray values from the 
reconstructed volume along the rays from a virtual source 

30 that intersect the modeled segment located between the 
virtual source and a virtual image plane. 
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In yet another aspect in accordance with principles of 
the present invention, the step of acquiring at least two x- 
ray images of the vascular structure is accomplished with a 
rotational acquisition using a C-arm x-ray system, the 
5 rotational acquisition at an angular velocity of at least 30° 
per second for a scan period of at least four seconds and 
less than six seconds. More preferably, the angular velocity 
is at least 55" during the scan period. 

An apparatus in accordance with principles of the 

10 present invention includes a support for an object affected 

by periodic motion to be imaged in an imaging region of an x- 
ray system. An x-ray system acquiries at least two x-ray 
images of the vascular structure. And ECG device obtains 
indicia of the phases of periodic motion affecting the 

15 object. A processor is included for correlating the indicia 
with each of the at least two x-ray images. A display is 
used for selecting at least two x-ray images from a similar 
phase of periodic motion- A processor is included to 
generate a reconstruction of a three dimensional modeled 

20 segment of a region of interest in the vascular structure, 

the modeled segment is reconstructed using the selected x-ray 
images from a similar phase of periodic motion, the region of 
interest only a portion of the imaged vascular structure. 
The processor generates a three dimensional voliometric 

25 reconstruction of a vascular structure larger than the 

modeled segment. In addition, the processor combines the 
modeled segment of interest and the volumetric reconstruction 
of the larger vascular structure- A display provides images 
of the generated reconstructions, the combined reconstructed 

30 . vascular model and volumetric reconstruction as well as other 
maps, models etc. in human readable form- 
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In a more limited aspect of the present invention, the 
processor generates a voliometric reconstruction of a first 
portion of the vascular structure reconstructed using x-ray 
images from all of the phases of periodic motion and 
5 generates a volumetric reconstruction for a second portion of 
the vascular structure with gated images from a similar phase 
of the periodic motion . 

In another aspect of an apparatus applying principles of 
the present invention the processor computes an overlap map 

10 of the vascular structure and modeled segment of interest, 
the overlap map computed by integrating all the gray values 
from the reconstructed volume along the rays from a virtual 
source that intersect the modeled segment located between the 
virtual source and a virtual image plane. 

15 The claimed method and apparatus combines the 

complementary features of volume and surface based model 
reconstruction techniques. The rotational acquisition is used 
to minimize contrast medium and X-ray exposure. The pro- 
jections that correspond to the same phase of the cardiac 

20 cycle are used to create an accurate surface-based model of 

the coronary segment of interest and to create an overview of 
the main coronary vessels. 

The presented novel method includes creation of a 
surface model of a selected segment of interest of the 
25 cardiac arterial tree based on projections from a rotational 
acquisition. The method has the potential to build a model 

based on every projection that captures the heart at the same 
phase. The combination of the partial model of the segment of 
interest with the automatic volumetric reconstruction of the 
30 entire coronary tree using all of the acquired projections, 
without regard to the cardiac phase, is used for 
visualization of the main vessels of the coronary tree and to 
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determine the ^optimal working view* with reduced vessel 
overlap and without the need to manually create a surface 
model of the entire coronary tree, thereby reducing 
reconstruction time making the process more useful for real 
5 time clinical application. 

Alternatively the method may use only the acquired 
projections that are around the same phase of the heart that 
are used in generating the model. This results in a displayed 
coronary tree structure that appears clearer but without the 
10 benefit of a representation of vessel position as a result of 
coronary motion. 

An apparatus and method applying principles of the 
present invention provides the foregoing and other features 
hereinafter described and particularly pointed out in the 

15 claims. The following description and accompanying drawings 
set forth certain illustrative embodiments applying 
principles of the present invention. It is to be appreciated 
that different embodiments applying principles of the 
invention may take form in various components and 

20 arrangements of components. These described embodiments 

being indicative of but a few of the various ways in which 
the principles of the invention may be employed. The 
drawings are only for the purpose of illustrating a preferred 
embodiment of an apparatus applying principles of the present 

25 invention and are not to be construed as limiting the 
invention. 

The foregoing and other features and advantages of the 
present invention will become apparent to those skilled in 
the art to which the present invention relates upon 
30 ' consideration of the following detailed description of a 
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preferred embodiment of the invention with reference to the 
accompanying drawings wherein: 

FIG. 1 is a schematic representation of an diagnostic 
imaging apparatus in accordance with principles of the 
5 present inventions- 
Fig. 2 is a flow diagram illustrating aspects of a 
method in accordance with principles of the present 
invention; 

Fig- 3 is a schematic representation of a process in 
10 accordance with principles of the present invention to 
generate an overlap map; and 

Fig. 4 is a flow diagram illustrating an alternate 
functional arrangement in accordance with principles of the 
present invention . 

15 With reference to Figure 1, a diagnostic imaging 

apparatus 20 includes an X-ray device 22, a video monitor 
bank 24 and diagnostic imaging control console 26. The x-ray 
device is a rotational angiography device comprising a C~arm 
system 30. A base 32 includes vertical column member 34 that 

20 supports the C-arm system 30. An arm 36 is movably attached 
to the vertical column member 43 for supporting a C-airm 38 
that exhibits an isocenter I. The C-arm 38 has an X-ray image 
pickup system formed by an X-ray source 4 0 and an X-ray 
receiver 42 that are respectively mounted at the opposite 

25 ends of the C-arm 38- The X-ray receiver 42 can be an image 
amplifier camera system or can also be a solid-state 
detector. The X-ray source 40 and the X-ray receiver 42 are 
arranged relative to one another such that a central beam of 
an X-rays emanating from the X-ray source 40 is incident 

30 approximately centrally on the X-ray receiver 42- 
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The C-arm 38 is motor-adjustable in the direction of the 
double arrow 44 along its circumference. The arm 36 can be 
rotated as shown by doioble arrow 4 6 by motor or manually 
around an axis A (angulation) - The C-arm is movable by motor 
5 or manually in a direction generally parallel with the axis A 
as shown by the double arrow 46. In addition, the C-arm 38 is 
vertically adjustable along the vertical column member 34 as 
shown by double arrow 48. Each of the moveable components are 
suitably monitored by known position encoders which provide 

10 the position of each component within the appropriate 

coordinate system for use by the imaging system control 
console 26. All of the motions described herein are capable 
of being carried out automatically by motor, or other means 
of causing the desired movement such as pneumatically, 

15 hydraulically etc., under the direction of the imaging system 
control console 26. 

The C-arm X-ray system 30 is provided for the generation 
of 3D images of a body area of an object 50 lying on a 
patient support (not shown) . A suitable x-ray device is a 

20 Philips Integris Allura® with a 12 inch monoplane. Philips 
Medical Systems an operating subsidiary of Koninklijke 
Philips Electronics N.V. having a principle place of business 
in Best, The Netherlands. 

Since periodic cardiac motion causes the coronary 

25 vessels in the arterial tree to move, a reconstructed 3-D 
image yields only a rough representation of the coronary 
arterial tree. In "order to improve the images used for 
generating the model used for the Quantitative Coronary 
Analysis, an ECG (electrocardiogram) device 52 is utilized to 

30 provide periodic coronary cycle information associated with 
each image acquired during the Rotational Acquisition. In 
addition the ECG information is also useful for generation of 
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a partial volumetric reconstruction, as described below- The 
" ECG is operatively connected with the system control console 
26 and has one or more body electrodes 54 for obtaining 
desired periodic physiological signals. 
5 The imaging system control console 26 coordinates the 

operation of the diagnostic imaging system 20. All of the 
control and imaging processing functions in the illustrated 
components and systems can be performed by known computer 
based systems having an operable complement of component 
10 systems such as suitable processors 60, memory and storage, 
input, output and data communications capabilities. A 
suitable workstation for the control console is an Octane 
system by Silicone Graphics Inc., with offices in Mountain 
View, California. 

15 An operator interface 90 includes input and output 

devices suitably in communication with the control console 2 6 
such as a keyboard 62, a touch screen monitor 64, a mouse 66> 
a joystick (not shown) , a track ball 68 as well as other 
input apparatus or devices to provide operator instructions 

20 to control the imaging system and generate the model of the 
arterial segment of interest in the vascular structure . 

Image processing and reconstruction circuitry in the 
control console 2 6 processes the output signals of the x-ray 
device 22 as it provides signals during an examination into 

25 an image representation. The image representation may be 
displayed on a video monitor, stored in computer memory, 
stored on tape or disk for later recall, further processed, 
or the like. 

Turning To Fig. 2, a method in accordance with the 
30 present invention is described for use with the apparatus 
described above. In step 100, the diagnostic imaging 
apparatus is initialized and calibrated according to known 
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methods. In step 102, the system parameters for the C-arm 
system such as positions of components and other conditions 
necessary for system operation are collected and provided to 
the system control 26. In step 104, the system acquires the 
5 images rotational acquisition and correlated ECG data for 
each acquired image. The C-arm is placed near the head- 
position of the table in order to perform a calibrated 
propeller acquisition from 120 RAO (right side of patient) to 
120 LAO (left side of patient) . The C-arm rotates with 55 
10 degrees per second and acquires images at 30 frames per 
second during 4 seconds while the patient holds his/her 

breath. A total of 8-12 cc's of contrast medium (Omnipaque 
350) is injected during the acquisition and started when the 

gantry starts rotating. This acquisition has been calibrated 

15 in step 100 and therefore, the individual projections can be 
corrected for the earth magnetic field and pincushion 
distortion and transformed to a common world coordinate 
system. The projections and the corresponding ECG 
information is stored in memory for selective recall and 

20 display on one of the monitors 24 or 64. 

Next, in step 106, a sub-set of the images acquired 
during the rotational acquisition are displayed on the 
monitor bank 24 or system monitor 64 for use by the 
clinician. In step 108, based on the ECG signal, the 

25 projections that correspond to the same selected phase of the 
cardiac cycle are chosen to model the segment of the artery 
of interest. Usually 2 projected images from different angles 
are selected. The modeled segment of interest generated from 
acquisitions from the same portion of the cardiac cycle 

30 provide the necessary accuracy to allow the clinician to mark 
the displayed image and have accurate measurements generated 



15 



wo 03/101300 



PCT/US03/17719 



from the model that can be used for accurate selection of the 
proper stent or other interventional procedure . 

It is to be appreciated that the entire coronary 
vascular structure is not, modeled- Modeling only the 
5 selected segment of interest, e.g. the stenosis or lesion 
which requires the improved accuracy for the Quantitative 
Coronary Analysis to be measured for the upcoming 
interventional procedure, reduces the computational time to 
construct the necessary maps described below such that the 
10 diagnostic process and interventional procedure may be 

carried out in real time with the patient in the surgical 
suite . 

Once the desired images are selected from the rotational 
acquisitions, the process moves to step 110 where the model 

15 segments are generated. A suitable method for generating the 
model segments is described in United States Patent No. 
6,501,848 entitled ^^Method and Apparatus for Three- 
dimensional Reconstruction of Coronary Vessels from 
Angiographic Images and Analytical Techniques Applied 

20 Thereto'' issued to Carroll et a,, which in incorporated in 
its entirety by reference herein. 

To construct the central axis of the modeled segment, 
the projections chosen above in step 108 include views in 
which the segment of interest is clearly visible and the 

25 angle between the two projections is around 90 degrees. Since 
on average five to six heartbeats are captured during the 
run, the user generally has a number of very desirable views 
from which to pick the corresponding projections with a clear 
view of the lesion. The central axis of the arterial segment 

30 is manually identified by the clinician by adding point pairs 
on the displayed image using any of the user interfaces 62, 
64, 66, 68 for the system controller 26. A point in 
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projection A yields an epipolar line in projection which 
is the ray coming from projection A towards the X-ray source. 
Once a point is set in one projection^ the user defines the 
corresponding point along the epipolar line in the other 
5 projection to construct the central axis of the model. By 

using this ^epipolar constraint', the corresponding epipolar 
lines intersect and define a 3D axis point. Once the axis of 
the lesion is created, the clinician delineates the borders 
of the lumen in every projection that corresponds to the same 

10 phase of the heart in order to create and refine a surface 

model of the liomen. The average time to create a model for a 
segment of interest, once two corresponding projections are 
selected, is less than a minute on a SGI Octane. 

Proceeding to step 112, the volumetric reconstruction is 

15 generated. An adapted version of the Feldkamp backpro j ection 
algorithm is applied to the acquired projections for the 
volumetric reconstruction. The projections are weighted 
according to the speed of time C-arm such that the 
projections acquired at a constant rotation speed have a 

20 higher weight than those acquired during startup and slowing 
down of the C-arm. If all available projections are used, the 
reconstruction contains information about several phases of 
the heart simultaneously, and is often difficult to 
interpret. If only those projections that correspond to the 

25 same time point in the cardiac cycle used to build the 

surface model, the reconstructed voliime is better suited for 
simultaneous visualization and inspection of the artery. 

The optimal view that is chosen by the clinician to 
perform the intervention is based on a combination of the 

30 amount of foreshortening of the segment of interest and the 
overlap of other coronary vascular structures. An optimal 
view map is so that only those C-arm gantry angles are 
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computed that can actually be achieved for use in the 
interventional phase of patient treatment. In general^ the 
cranial (toward the head of the patient) and caudal (toward 
the feet of the patient) angulation should not exceed 30 
5 degrees from the anterior-posterior (AP) plane since the 

image intensifier contact the patient otherwise. The RAO and 
LAO thresholds can be set according to the capabilities of 
the imaging system and are currently set, for the Philips 
system described above, to a maximum of 60 degrees from the 

10 AP plane to avoid collision of the image intensifier with the 
operating table or the patient. The different values within 
the predefined range of angles are visualized using an 
interactive color coded map. The different portions of such a 
map correspond to different positions e.g., quadrant or area 

15 of the map for a RAO position another corresponds to a 

cranial position (bottom with caudal) . The clinician has the 
opportunity to inspect the different viewing angles in real 
time by clicking at a position in the map displayed on the 
monitor 64 and to toggle between the different maps. The 

20 amount of foreshortening and overlap are given in the lower 
right hand side of the display. 

In step 114, the foreshortening map is computed by 
comparing the length of the modeled segment to the projected 
length of this segment as if it were viewed from a typical 

25 viewpoint as defined by the range of angles in the map. 

Next, in step 116, a novel method of computing the 
overlap map is carried out by the control console 26. 
Although the volumetric reconstruction of all the acquired 
projections yields only a rough representation of the 

30 coronaries, it does provide useful information such as the 
location of the vessels and other objects (e.g. the spine. 
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ribs and pacemaker or ECG-leads) during all the phases of the 
heart . 

Turning now to Fig. 3, a simplified example is 
schematically illustrated of the novel method, in accordance 
5 with the present invention, for determining the overlap map. 
The overlap map is computed by taking the integral of all the 
gray values from the reconstructed volume along the rays that 
intersect the modeled segment- A virtual x-ray source 140 has 
a plurality of rays 142 corresponding to a volumetric 

10 reconstruction 144. A lumen 146 is represented in the 

reconstruction with a modeled segment 148 of interest such as 
a stenosis or lesion. An overlap area 149 is where another 
portion of the lumen or another structure intercepts a ray 
that has also passed through the modeled segment 148. A 

15 planar representation of the sum of gray values 150 has a 

corresponding lumen representation 14 6a of the lumen 14 6 and 
its model segment 148 shown as 14 8a in representation 150. 
An illustration of the sum of gray values on the model 
segment 152 is shown to include an area 154 defining the^ 

20 overlap area of the modeled segment and the lumen. The four 
dashed lines 156a,b,c,d schematically illustrate a 
visualization sxxm of gray values along a ray that are 
projected on the model. This represents the concept of the 
overlap of the vascular structure creating a distribution 

25 which provides the representation in the graph 158 with 
intensity illustrated on the axis 160. The graph 158 
illustrates the intensities (=sxjm of gray values along the 
ray) over a single line. The final value is the integral or 
sum of all the sum of gray values in the projection of the 

30 modeled segment.. This gray value integration technique is 
efficiently implemented in OpenGL using space-leaping 
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techniques and the Stencil Buffer in coiabination with the 
special purpose graphics hardware of the SGI Octane - 
Next^ in step 118, the optimal view map can he 
constructed by taking a weighted sum of normalized values of 
5 the foreshortening map and overlap map. For a practical 

implementation and efficient computation, the vessel overlap 
is generally estimated for gantry angles where the vessel 
foreshortening has been reduced or minimized in order to 
obtain the resultant optimal working view. 

10 In step 120, the combined reconstructed image including 

the voliimetric reconstruction and modeled segment, the 
optimal view map and system parameters to locate the gantry 
for a selected optimiom view are displayed on monitors in the 
monitor bank 24 or control console 26. The process 

15 terminates in step 122. 

Turning to Fig. 4, an alternate method applying 
principles of the present invention is illustrated. In step 
200, the diagnostic imaging apparatus is initialized and 
calibrated according to known methods. In step 202, the 

20 system parameters for the C-arm system such as positions of 
components and other conditions necessary for system 
operation are collected and provided to the system control 
26. In step 204, the system acquires the images rotational 
acquisition and correlated ECG data for each acquired image. 

25 The C-arm is placed near the head-position of the table in 

order to perform a calibrated propeller acquisition from 120 
RAO (right side of patient) to 120 LAO (left side of 
patient) . The C-arm rotates with 55 degrees per second and 
acquires images at 30 frames per second during 4 seconds 

30 while the patient holds his/her breath. A total of 8-12 cc's 
. of contrast medium (Omnipaque 350) is injected during the 
acquisition and started when the gantry starts rotating. This 
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acquisition has been calibrated in step 100 and therefore, 
the individual projections can be corrected for the earth 
magnetic field and pincushion distortion and transformed to a 
common world coordinate system- The projections and the 
5 corresponding ECG information is stored in memory for 

selective recall and display on one of the monitors 24 or 64. 

In step 205, the system is adapted to generate both or 
either of (i) a full volumetric reconstruction in step 212, 
similar to step 112 above in Fig. 2 and (ii) a cardiac gated 

10 volumetric reconstruction of a portion of the arterial tree, 
for example, including the main branches for the tree which 
includes the modeled segment of interest. In step 205, if 
the system is set to generate only the full volumetric 
reconstruction in step 112, the process is similar to that in 

15 Fig. 2. Alternatively, the system may only generate a 

partial volumetric reconstruction using only images from the 
same portion of the cardiac cycle, ECG gained images, to 
generate a more clear partial cardiac arterial tree which 
includes the modeled segment of interest. In addition, the 

20 system may generate both a full volumetric reconstruction 
using all acquired images which will provide information 
illustrating aspects of the cardiac motion and generate the 
partial volumetric reconstruction which includes the modeled 
segment . 

25 Next, in step 208, a sub-set of the images acquired 

during the rotational acquisition are displayed on the 
monitor bank 24 or system monitor 64 for use by the 
clinician. Based on the ECG signal, the projections that 
correspond to the same selected phase of the cardiac cycle 

30 are chosen to model the segment of the artery of interest as 
well as provide- the images for the partial reconstruction 
that is generated in step 207. Usually 2 projected images 



21 



wo 03/101300 



PCT/US03/17719 



from different angles are selected for the model and all of 
the images from the same portion of the cardiac cycle, 4-6 
gated images in a typical scan in this method, are used for 
the partial voliametric reconstruction. The modeled segment 
5 of interest generated from acquisitions from the same portion 
of the cardiac cycle provide the necessary accuracy to allow 
the clinician to mark the displayed image and have accurate 
■measurements generated from the model that can be used for 
accurate selection of the proper stent or other 

10 interventional procedure. 

As discussed above for Fig. 2, it is to be appreciated 
that the entire coronary vascular structure is not modeled. 
Modeling only the selected segment of interest^ ©-9- the 
stenosis or lesion which requires the improved accuracy for 

15 the Quantitative Coronary Analysis to be measured for the 

upcoming interventional procedure, reduces the computational 
time to construct the necessary maps described below such 
that the diagnostic process and interventional procedure may 
be carried out in real time with the patient in the surgical 

20 suite. 

Once the desired images are selected from the rotational 
acquisitions, the process moves to step 210 where the model 
segments are generated. The process in step 210 is similar 
to that described above for step 110 in Fig. 2. 

25 As mentioned above, in step 112, the volumetric 

reconstruction is generated. An adapted version of the 
Feldkamp backproj ection algorithm is applied to all of the 
acquired projections for the volumetric reconstruction. The 
projections are weighted according to the speed of time C-arm 

30 such that the projections acquired at a constant rotation 

speed have a higher weight than those acquired during startup 
and slowing down of the C-arm. Since all available 
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projections are used, the reconstruction contains information 
about several phases of the heart simultaneously, and is 
often difficult to interpret. Also described above, in step 
207, the Feldkamp Backprojection algorithm is applied to only 
5 . those projections that correspond to the same time point in 
the cardiac cycle used to build the surface model, the 
partial reconstructed volume is better suited for 
simultaneous visualization and inspection of the artery with 
the modeled segment of interest- The volumetric 

10 reconstruction information from steps 207 and/or 212 is 
provided to step 216. 

In step 216, the overlap map is computed as described 
above for step 116 of Fig. 2 and is carried out by the 
control console 26. Although the volumetric reconstruction 

15 of all the acquired projections yields only a rough 

representation of the coronaries, it does provide useful 
information such as the location of the vessels and other 
objects (e.g. the spine, ribs and pacemaker or ECG-leads) 
during all the phases of the heart. If the partial 

2 0 reconstruction is included the branch of the coronary 
arterial tree is more clearly illustrated. 

In step 214, the foreshortening map is computed by 
comparing the length of the modeled segment to the projected 
length of this segment as if it were viewed from a typical 

25 viewpoint as defined by the range of angles in the map. 

Next, in step 218, the optimal view map can he 
constructed by taking a weighted sum of normalized values of 
the foreshortening map and overlap map. For a practical 
implementation and efficient computation, the vessel overlap 

30 is generally estimated for gantry angles where the vessel 
foreshortening has been reduced or minimized in order to 
obtain the resultant optimal working view. 
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In step 220^ the system displays the combined 
reconstructed image including one or both of the full and/or 
partial volumetric reconstruction and modeled segment^ the 
optimal view map and system parameters to locate the gantry 
5 for a selected optimum view. The images models and data are 
displayed on monitors in the monitor bank 24 or the control 
console 26. The process terminates in step 222. 

A hybrid reconstruction scheme including principles of 
the present invention is described that provides the 

10 clinician with three-dimensional quantitative measurement of 
an accurately modeled arterial segment and an optimal working 
view throughout the cardiac cycle. This information can be 
provided in less than a minute and used to plan the 
intervention while the patient is still on the 

15 catheterization table. 

While a particular feature of the invention may have 
been described above with respect to only one of the 
illustrated embodiments, such features may be combined with 
one or more other features of other embodiments, as may be 

20 desired and advantageous for any given particular 
application. 

From the above description of the invention, those 
skilled in the art will perceive improvements, changes and 
modification. For example, other projections captured at the 

25 same phase can be used to refine the surface model of the 

segment of interest. The projections that are approximately 
of the same phase may be used to visualize the coronary tree, 
while, all the acquired projections can be used to compute 
the overlap map. It is to be appreciated that either the 

30 full or partial volumetric reconstruction can be used to 

compute the overlap map or be used to generate images on one 
of the monitors for visualization of the hybrid coronary 3-D 
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structure. Such improvements, changes and modification 
within the skill of the art are intended to be covered by the 
appended claims . 
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CLAIMS : 

1. A method of generating a hybrid three dimensional 
reconstruction of a vascular structure affected by periodic 
motion,, the method comprising: 
5 placing an object affected by periodic motion to be 

imaged in an imaging region of an x-ray system, the 
object having a vascular structure; 

acquiring at least two x-ray images of the vascular 
structure; 

10 obtaining indicia of the phases of periodic motion 

and correlating the indicia with each of the at least 
two x-ray images; 

selecting at least two x-ray images from a similar 
phase of periodic motion; 

15 generating a three dimensional modeled segment of a 

region of interest in the vascular structure, the 
modeled segment reconstructed using the selected x-ray 
images from a similar phase of periodic motion, the - 
region of interest only a portion of the imaged vascular 

20 structure; 

generating a three dimensional volumetric 
reconstruction of a vascular structure larger than the 
modeled segment; 

combining the modeled segment of interest and the 

25 volumetric reconstruction of the larger vascular 

structure; and 

displaying in human readable form the combined 
reconstructed vascular model and volumetric 
reconstruction . 

30 
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2 . The method of claim 1 wherein the volumetric 
reconstruction includes using x-ray images from all of the 
phases of the periodic motion. 

5 3 . The method of claim 1 wherein the volumetric 

reconstruction is generated with gated images from a similar 
phase of periodic motion. 

4 . The method of claim 2 wherein the volumetric 
10 reconstruction includes a first portion of the vascular 

structure reconstructed using x-ray images from all of the 
phases of periodic motion and a second portion of the 
vascular structure reconstructed with gated images from a 
similar phase of the periodic motion. 

15 

5 . The method of claim 4 wherein the reconstruction of both 
of the first portion and second portion of the vascular 
structure are combined with the modeled segment of interest. 

20 6. The method of claim 1 wherein the obtained indicia of 

phases of periodic motion is representative of cardiac motion 
and is provided by an ECG signal. 

7 . The method of claim 1 including computation of an 

25 overlap map, the overlap map computed by integrating all the 
gray values from the reconstructed volume along the rays from 
a virtual source that intersect the modeled segment located 
between the virtual source and a virtual image plane. 

8 . The method of claim 1 wherein acquiring at least two x- 
30 ray images of the vascular structure is accomplished with a 

rotational acquisition using a C-arm x-ray system, the 
rotational acquisition at an angular velocity of at least 30° 
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per second for a scan period of at least four seconds and 
less than six seconds . 

9. The method of claim 8 wherein the angular velocity is at 
5 least 55° during the scan period. 

10. An apparatus for generating a hybrid three dimensional 
reconstruction of a vascular structure affected by periodic 
motion^ the apparatus comprising: 

10 means for supporting an object affected by periodic 

motion to be imaged in an imaging region of an x-ray 
system, the object having a vascular structure; 

means for acquiring at least two x-ray images of 
the vascular structure; 

15 means for obtaining indicia of the phases of 

periodic motion and correlating the indicia with each of 
the at least two x-ray images; 

means for selecting at least two x-ray images from 
a similar phase of periodic motion; 

20 means for generating a three dimensional modeled 

segment of a region of interest in the vascular 
structure, the modeled segment reconstructed using the 
selected x-ray images from a similar phase of periodic 
motion, the region of interest only a portion of the 

25 imaged vascular structure; 

means for generating a three dimensional volumetric 
reconstruction of a vascular structure larger than the 
modeled segment; 

means for combining the modeled segment of interest 

30 and the volumetric reconstruction of the larger vascular 

structure; and 
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means for displaying in hioman readable form the 
combined reconstructed vascular model and volumetric 
reconstruction. 

5 11. The apparatus of claim 10 wherein the means for 

volumetric reconstruction uses x-ray images from all of the 
phases of the periodic motion. 

12. The apparatus of claim 10 wherein the means for 
10' voliometric reconstruction generates the volumetric 

reconstruction with gated images from a similar phase of 
periodic motion. 

13. The apparatus of claim 11 wherein the means for 

15 volumetric reconstruction includes means for generating a 

volumetric reconstruction of a first portion of the vascular 
structure reconstructed using x-ray images from all of the 
phases of periodic motion and means for generating a 
volumetric reconstruction for a second portion of the 

20 vascular structure with gated images from a similar phase of 
the periodic motion. 

14. The apparatus of claim 13 wherein the means for 
combining the modeled segment of interest with the volumetric 

25 reconstruction includes means for combining both of the first 
portion and second portion of the vascular structure with the 
modeled segment of interest. 

15. The apparatus of claim 10 including an ECG monitor to 
30 obtain the indicia of phases of periodic motion . 
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16. The apparatus of claim 1 including means for computation 
of an overlap map,, the overlap map computed by integrating 
all the gray values from the reconstructed volume along the 
rays from a virtual source that intersect the modeled segment 

5 located between the virtual source and a virtual image plane. 

17. A method of generating a hybrid three dimensional 
reconstruction of a vascular structure affected by periodic 
motion;, the method comprising: 

10 placing an object affected by periodic motion to be 

imaged in an imaging region of an x-ray system, the 

object having a vascular structure; 

acquiring a rotational acquisition of x-ray images 

of the vascular structure; 
15 obtaining indicia of the phases of periodic motion 

and correlating the indicia with each of the at least 

two x-ray images; 

selecting at least two x-ray images from a similar 

phase of periodic motion; 
20 generating a three dimensional modeled segment of a 

region of interest in the vascular structure, the 

modeled segment reconstructed using the . selected x-ray 

images from a similar phase of periodic motion, the 

region of interest only a portion of the imaged vascular 
25 structure; 

generating a three dimensional volumetric 

reconstruction of a vascular structure larger than the 

modeled segment; 

combining the modeled segment of interest and the 
30 volumetric reconstruction of the larger vascular 

structure; and 
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displaying in hioman readable form the combined 
reconstructed vascular model and volumetric 
reconstruction . 

5 18. The method of claim 17 wherein the rotational 

acquisition acquires images at an angular velocity of at 
least 30° per second for a scan period of at least four 
seconds and less than six seconds - 

10 19. The method of claim 8 wherein the angular velocity is at 
least 55° during the scan period. 

20. The method of claim 17 wherein the volumetric 
reconstruction includes using x-ray images from all of the 

15 phases of the periodic motion. 

21. The method of claim 17 wherein the volumetric 
reconstruction is generated with gated images from a similar 
phase of periodic motion. 

20 

22. The method of claim 17 wherein the volumetric 
reconstruction includes a first portion of the vascular 
structure reconstructed using x-ray images from all of the 
phases of periodic motion and a second portion of the 

25 vascular structure reconstructed with gated images from a 
similar phase of the periodic motion. 

23. The method of claim 22 wherein the reconstruction of 
both of the first portion and second portion of the vascular 

30 structure are combined with the modeled segment of interest. 
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24. The method of claim 17 wherein the obtained indicia of 
phases of periodic motion is representative of cardiac motion 
and is provided by an ECG signal. 

5 25. The method of claim 17 including computation of an 

overlap map, the overlap map computed by integrating all the 
gray values from the reconstructed volume along the rays from 
a virtual source that intersect the modeled segment located 
between the virtual source and a virtual image plane. 
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